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Abstract. 

We consider correlation functions for the Wess-Zumino-Witten model on the torus 
with the insertion of a Cartan element; mathematically this means that we consider 
the function of the form F = Tr('I>i(2;i) . . . ^n{zn)q^^e^) where <I>i are intertwiners 
between Verma modules and evaluation modules over an afRne Lie algebra g, d is 
the grading operator in a Verma module and h is in the Cartan subalgebra of g. We 
derive a system of differential equations satisfied by such a function. In particular, 
the calculation oi q-^F yields a parabolic second order PDE closely related to the 
heat equation on the compact Lie group corresponding to g (cf. [Ber]). We consider 
in detail the case n = 1, g = s[2. In this case we get the following differential 
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equation (g = e'"'^): (^-27ri(/C + 2) ^ + f^j F = {m{m + l)p{x + ^) + c)F , which 

for K = —2 (critical level) becomes Lame equation. For the case m £ Z we derive 
integral formulas for F and find their asymptotics as K — > —2, thus recovering 
classical Lame functions. 



Introduction. 

We start with consideration of the Wess-Zumino-Witten model of conformal field 
theory on a torus. This is, we consider an affine Lie algebra g corresponding to some 
simple finite-dimensional Lie algebra g. For technical reasons, it is more convenient 
to work with a twisted realization of g. Next, we consider Verma modules Mx^k over 
g. If y is a representation of the finite-dimensional algebra g then by definition a 
vertex operator $(z): M\^k — > M^^k is an operator valued formal Laurent series 
in z satisfying the following commutation relations with the elements of g: 

$(z)a ® t"" = ((a ® t"^) ® 1 + z""! ® a) 

Let Mx-^k, i = . . .n be a collection of Verma modules such that Aq = A^, and 
$^(^1) : Mxi^k Mxi_j^^k ® Vi be vertex operators. Then we can consider the 
following "correlation function on the torus" : 



where d is the grading operator in Verma modules^ and h e f)K. This function 
takes values in the module F = Fi (8) . . . (8) and it is the main object of our study. 

Our first goal is to derive differential equations for JF. We compute -^J-' using 
the same technique as for the usual Knizhnik-Zamolodchikov equations (see [TK], 
[FR]). However, this system of equations (Theorem 3.1) is not closed: it has the 
form 

d d 
Zi-Q^T =Ai{zi... Zn)^ + J2 ^*^^')^* 

where Ai arc some operators in V and the sum is taken over an orthonormal basis 
xi in f). Since we do not have any information about this system does not 

allow us to determine JF. This system of equations in another form appeared first 
in the papers of Bernard ( [Ber] ) . 

However, we can get additional information studying dependence on the modular 
parameter q. This dependence takes a simpler form if we renormalize JF, considering 
J-'/J-'oj where J-'o = Trl^o k {q~^^'^)- For this function we have the following equation 
(Theorem 4.1): 

where B is some operator in V depending on Zj, q. 

This equation is especially interesting for n = and n = 1, since in these cases 
there is no dependence on z. For n = we obtain that T jT^ satisfies the heat 
equation on the maximal torus T d G. This result was first obtained by Bernard; 
however, he worked in the untwisted setting which gave him heat equation on the 
group rather than on the torus. 

We concentrate our attention on the case n = 1, g = and take V to be 2m + 1- 
dimensional irreducible module over sl2- In this case the intertwiner (^{z): M\^k 
Mx^k ® V exists and is unique up to a constant for generic values of k, A. Therefore 
we can apply the results of the previous sections. It will be convenient to identify 
f)R ~ R and consider the following function of one real variable x: 

_ TV|m,,, ($(^)Q-^e-^"^) 
where h = ( }. )Gs[2. Then the results of Section 4 show that thus defined 



^0 -1^ 

F satisfies the following partial differential equation (eq. (5.1)): 

-47ri(A; + 1)^ + ^ J F = m(m + 1) + ^) + c) F, 



where q = 



-"^We use the symbol d for the grading operator in twisted reahzation, reserving the standard 
notation d for the untwisted grading operator, see Section 1. 
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This equation has many nice properties and is worth studying in itself; it is a 
non-stationary Schrodinger equation. Equations of this type were considered and 
treated in [Kri] and references therein. What is particularly interesting is that in 
the limit /c — > — 1 this equation becomes the Lame equation with precisely the 
same expression of the coefficient to p as in the classical form (i.e., m{m + 1), cf. 
[WW]). Note that the case = — 1 in the twisted setting corresponds to K = — 2 
in the untwisted one, i.e. to the critical level. Therefore we should expect that the 
asymptotics of the correlation functions on the torus at the critical level limit are 
closely related to the solutions of Lame equation. 

We make this relation more explicit. First we derive integral formulas for the 
function F defined above using the Wakimoto realization of Verma modules over 
5(2- The answer is given by Theorem 5.1 and has the following form: 

F{x,q) = J l{C\T)'/^<l>{x,C\r)dC,...dCm 

A 

where both and (f){x, Q are expressed in terms of theta functions and exponents 

and K = 2{k + 1). 

Now we can find the asymptotics of these functions as k — > 0. To do it we need 
the saddle-point method which says that under certain assumptions this integral 
has the following asymptotics: 

F{x,q) ^ C{T)K^/H{Co\ry/''<j>{x,Co\r), 

where (0 is a critical point for the function 1{C). Substituting this asymptotics in 
equation (5.1) we finally see that the function 

f{x) = (pixXo) 

is a solution of the Lame equation. Note that we cannot find the critical point 
Co explicitly; we can only write a system of equations defining it; this system of 
equations will also depend on the highest weight A of the Verma module. 

It turns out that in this way we can get a generic eigenfunction of the Lame 
operator; in particular, we can get all the doubly periodic eigenfunctions which 
are obtained for integer A. For example, for m = 1 we show how to get the 
classical solutions f{x) — dn 2Ka;,sn 2Ka;,cn 2'Kx, which correspond to even (the 
first solution) and odd (two latter solutions) values of A. 

However, it is not very easy to check the applicability of the saddle-point method. 
For the case m = 1, we can check it in one special case (g G M, A is small) and then 
use analytic continuation arguments to show that in fact every critical point gives 
rise to an eigenfunction of Lame equation. In the general case one could repeat 
similar arguments to prove that again any critical point gives a solution of Lame 
eqution; however, we do not follow this way; instead, we say that it can be checked 
by direct calculation which is in fact carried out in the book by Whittaker and 
Watson [WW]. 

Remark. The asymptotics of solutions of the KZ equations at the critical level were 
computed by means of the saddle-point method by Reshetikhin and Varchenko[RV]. 
They found that the leading term of these asymptotics is a solution of the Bethe 
equations well known in physics. 
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1. A twisted realization of afflne Lie algebras. 

Let g be a finite dimensional simple Lie algebra over C of rank r. Denote by 
<, > the standard invariant form on q with respect to which the longest root has 
length \/2. 

Let f) denote a Cartan subalgebra of Q. The form <, > defines a natural identi- 
fication X hx ioT X e i)* . We will use the notation <, > for the inner 
product in both t) and [)*. 

Let A+ be the set of positive roots of q. For a G A+, let G fla, /a €i 0_q, 
be such that < Ca, fa >— 1- Then [ca, fa] — ha- Also, let xj, 1 < i < r, be 
an orthonormal basis of with respect to the standard invariant form. Then the 
elements Cq,, Xi form a basis of g. 

Let g be the dual Coxeter number of g. Let p = ^ 'I2aeA+ ^- Then for a G 
A"*", < p, a >= \oi\ is the number of summands in the decomposition of a in the 
sum of simple positive roots. 

Let 7 = Ad e 9 be an inner automorphism of g. This automorphism is of 
order g. 

The action of 7 on root vectors is as follows: 7(60,) = e'^^'eQ, , 7(/a) = e~'"'/a 1 
where e = e^'^^l^ is a primitive ^-th root of unity. Note that since < |q!| < ^ — 1, 
7(a) = a if and only if a G I). 

Let = 0® C[t,t~^] © Cc be the affine Lie algebra associated with g. The 
commutation relations in this algebra are 

(1.1) [a(t) + Ac, hit) + /ic] = [a(t), h{t)\ + — { < a'{t)b{t) > r'^dt ■ c 

Stti J\t\=i 

for any two g-valued Laurent polynomials o(t), and complex numbers A, 

Define the twisted affine algebra subalgebra of consisting of all expres- 

sions a{t) + Ac with the property a{et) = 7(a(t)). The basis in this subalgebra is 
given by elements Ca ® /a ® Xi®^^^, c, for m G Z, a G A+. For 
brevity we write x[n\ for x <S)t"' and x for x ® t^. 

The representation theory of is quite parallel to that in the untwisted case. 
First of all, define the polarization of g^: 0^ = ® 0~ ® f) ® Cc. Here 0+ is the set 
of polynomials a{t) vanishing at 0, and 0~ is the set of polynomials a{t) vanishing 
at infinity. 

Next, define Verma modules over 0^. This is done exactly in the same way as 
for the untwisted affine algebra. Let A G f)* be a weight, and let be a complex 

^Note that throughout the paper we denote the complex number i — \f—\ by a roman "i", to 
distinguish it from the subscript which is itaUc. 
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number. Define Xx^k to be a one dimensional module over 0+ © {) ® Cc spanned by 
a vector v such that 0+ annihilates v, and cv = kv, hv = X{h)v, h & 1). Define the 
Verma module 

(1-2) ^*.'=K:e,^cc^>* 

Now define evaluation representations. Let be a module over g; we always 
assume that V has a weight decomposition such that weight subspaces are finite- 

27rihp 

dimensional. Define the operator C G End(F) by C = e s ; then Caw — j{a)Cw 
for any w & V , a & Q, and Cwq = wq, if wq is from the zero-weight subspace V{0). 

Let V{z) denote the space of V- valued Laurent polynomials in z, and let Vc{z) 
be the space of those polynomials which satisfy the equi variance condition w{ez) = 
Cw{z). 

The natural (pointwise) action of g on V{z) restricts to an action of on Vc{z). 
Note that the twisted affine algebra is isomorphic to untwisted one. More pre- 
cisely, we have the following lemma. 

Lemma 1.1 (see [PS], p. 36). The Lie algebras and g are isomorphic. Under this 
isomorphism, the Verma module M^,/; over is identified with the Verma module 
■M\+kp,gk over g, and evaluation representation Vc{z) over g^y is identified with the 
evaluation representation V{z) over g. 

Therefore, all the results about the representations of g-y can be as well ob- 
tained from the representations theory of g. However, use of the twisted algebra is 
technically more convenient to us. 

We can also extend the Lie algebra g^ adding to it the degree operator d which 
commutes with elements of g^ by [d,a{t)] = ta'{t),[d,c] = 0. We denote this 
extended Lie algebra by g^. Note that under the isomorphism of Lemma 1.1, 
d ^ gd+ hp, where d is the standard grading operator for g. 

Since every Verma module has a natural Z-gradation, we can uniquely extend 
the action of g^ in Mx^k to the action of g^ as soon as we define the action of d 
on the highest weight vector. Let us define it by dvx,k — — 2{k+i) '^A.fe- Again, this 
agrees with the untwisted case: if we identify Mx,k with the Verma module A4a,k 
over g, where A = A + kp, K = gk, then d is identified with gd + hp — '^2(k+i)' i 
where the action of d on Verma module M.x,k is defined so that on the highest 
level, (i = — ^2(K+g^^ • This gives us the following twisted analogue of Sugawara 
construction for d. 

Proposition 1.2 (see [E]). In every Verma module over g^, 
(1.3) 



d = 



~k~j:i I • + ^9\fa[-\o.\ -mg]:^-^^: Xj[mg\xj[-mg\ : j , 

m€Z yQ:eA+ j=l J 



k + 

where :: is the standard normal ordering 

■ ea[n\fa[-n\ := 



ea[n]fa[~n] n<0 
U[-n]ea[n] n>0 

(1.4) : h[n]h[-n] := \ , V^,. 7 n ' ^ ^ 

h[—n\h[n\ n> 
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We can also define the action of d in the evaluation representation Vc{z) by 
d = z-^. Thus Vc{z) becomes a g^-module. Again, under the isomorphism of the 
Lemma 1.1, d is identified with the operator gd+hp in the evaluation representation 
V{z) over g, where we define the action of d in V{z) hj d = z-j^.^ 

Let us introduce the twisted version of currents. Set 

Je„(^)= 5]e«[|a|+m^]-^-l"l— 
Mz)=J2h[mg]-z-^^-\ hei). 

mel. 

Thus by linearity we have defined Ja{z) for any a & g. 
Define the polarization of currents: 

m<0 

Jliz) = U[-\a\ + mg] ■ z\-\-^^-\ 

m<0 

(1.6) J+{z) ^^h-z-'+Y^ h[mg] ■ z-^^-\ /I e (). 

m<0 

This defines Jai^) ^or all a e q. Now set 

(1.7) J:(z) = J+{z)-Ja(z). 

Note that this polarization is not quite the same as the standard polarization 
of currents for the untwisted q, i.e. the isomorphism between q and does not 
match up these two polarizations. 

One can easily write commutation relations between currents; however, we won't 
need them. One thing we will need is the commutation relations with d. Namely, 
one can easily see that 

(1.8) Jtiz)q-^ = q-\-^J^iq-'z) 



(1.5) 



2. Twisted intertwiners and correlation functions on a torus 

We will be interested in intertwining operators ^{z) : Mx,k ~^ Mi,^k<^z^Vc{z), 
where the highest weight of F is /i, ® denotes the completed tensor product, and 
A is a complex number. 

Let zq be a nonzero complex number. Evaluation of the operator $(^) at the 
point Zq yields an operator ^(zq) : M\^k — A^i^,fe ® V-, where M denotes the 
completion of M with respect to the grading. 

^Note, however, that usually in the conformal field theory the action of d on the evaluation 
reperesentation is defined with some shift, namely: d — z-^ + 2{K-\-g) ' ^^ere C is the Casimir 
operator. We will note use this convention here. 
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From now on the notation ^(z) will mean the operator $ evaluated at the point 
z G C*. This will give us an opportunity to regard the operator $(-2) as an analytic 
function of z. This analytic function will be multivalued: ^{z) — z^^^{z), where 
$° is a single- valued function in C*, and A = '^"'"■^k+i)'^^ ■ 

Let u belong to the restricted dual module V*. Introduce the notation 
^u{z) = {l®u){^{z)). ^u{z) is an operator: Mx,k ^ M,,^k- 

The intertwining property for ^{z) can be written in the form 

(2.1) [a[m],^u{z)] = z'^^au{z). 

It is convenient to write the intertwining relation in terms of currents. 
Lemma 2.2. 

[Ji(C), ^u{z)] = r^x«(C/^)$e„.W, a e A+; 

(2.2) [jl{0,^^{z)]=C'x--{C/z)^f^^{z), ae A+, 



where 



At) = iii^. 

t9 



x"{t)=t-\ 
(2.3) x-''{t)=t 



1^ 



The identities marked with + make sense if \z\ > and those marked with — 
make sense if \z\ < 

Now we are ready to write down the twisted version of the operator Knizhnik- 
Zamolodchikov (KZ) equations. 

Theorem 2.3 (see[E]). The operator function ^u{z) satisfies the differential equa- 
tion 

{k + l)^$n(^) = Yl (•^et(^)^/«u(^) - '^f.u{z)J-Jz)) + 

(2.4) 

r 

E (4J^)*e„.(^) - ^e^u{z)JjJz)) + Y.{J+{Z)^.M^) - ^.,u{z)J-{z)). 

aeA+ 3=1 

In the paper [E] this proposition is proved in the case when y is a highest-weight 
module. However, one can check that the result is valid for any module with the 
weight decomposition. 

Remark. Note that since we have the isomorphism of Lemma 1.1, we can identify 
the intertwiners for with those for g. Namely, if $(C): -^a,k Mn,k <8) V]u is an 
intertwiner for 9, then the results of Section 1 show that 
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(2.5) 



$(2) =^-^^$(^9) 



will be an intertwiner for g^. 

The main object of our study will be the following twisted version of correlation 
functions on a torus. Let Vi, Vn be representations of Q, and let ^'^{zi) : M\^^k 
M\^_-^,k Fi, 1 < i < n be intertwining operators for (sometimes for brevity we 
write instead of $^.(2^)). If Aq = A„ = A, then we can define the following 
correlation function on a torus: 

(2.6) Tu„... ,uS^u ...,zn;h)= IV|m,.,«(^i) • . .KAzn)q-''e% 

where q G C^, \q\ < 1, h E i)R is an element from the compact form of f) (i.e., 
h{a) G iM for all ct G A). This trace is a formal power series in 2;i . . . Zn- However, 
it can be shown that it converges to an analytic function of Zi in the region l-^ij > 
. . . > \zn\ > \qzi\. Therefore, we will consider as a function of q^zi, . . . ,Zn,h 
with values in Vi . . . ^ or, equivalently, as a function of q, Zi with values in 
Vi ^ . . . ^ Vn ® C°^{i)-R). Our first goal is finding differential equations for this 
function. 

3. Differential equations for correlation functions. 

In this section we deduce the differential equations for the correlation function 
on the torus defined by (2.6). The basic tool will be the application of operator 
KZ equation (2.3), which implies: 

d 

[k + l)Zi-^^ J^ Ul,... ,Un i^lj • • • ) ^n) = 

+ J2 ^iit(^^..(j;j^o<.,(^.)-$u(^0"^/:(^o)---*>-^e'^) 

r 

(3.1) + J]^.TV($^ . . . {J+{z,)K,uM) - K,uM)J^.i^^)) ■ ■■^''Q-'e^). 
1=1 

Let us consider the first summand in (3.1). Using commutation relations of 
currents with intertwiners (2.2), we can move J^^ (zi) to the utmost left, and (zi) 
to the utmost right, which gives us the following expression for the first term in 
(3.1): 

(3.2) 

aeA+ ^ jy^i 

- z,Tr{Jl (z,)^' . . . $)^,^ (z,) . . . $-g-^e^) 

+z,Tr{^' . . . (z,) . . . J- {z,)q-^e^)^ . 
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Now, using the cyclic property of trace and the relation q ^e^jf^ {z) = qe'^^°'^ Jf^ {(iz)q 
we can rewrite (3.2) as follows: 



aeA+ ^ i/i 

- e'^(«)g^,Tr($i . . . {z^) . . . ^'^ J+ {qzi)q-^e^) 

+ e-'^(«)?-i^,IV( J- (q-'zi)^' . . . {Zi) . . . ^^q-^e^)^ 



0, we finally get the following expression for the first summand in (3.1) 



Repeating this procedure and taking into account that lim g z Jg {q "^z) 

m — >+cxD ™ 



(3-3) -EE ^"(?"^^^A.)e"^'^^"V,®7r,(^®e«)$, 

Q:eA+ mez,j 

where means that the sum is taken over all m G Z, j = 1 . . . n except m = 0, j = 
i. (We use the convention tt^ ® TVi{a ^ b) = 7ri{ba).) It is easy to see that the series 
fa,h — ^ x'^{q'^Zi/ Zj)e'^^^°'^ converges; it can be expressed via theta-functions, 

though we won't need it. 

Similarly, one shows that the second term in (3.1) equals 



(3-4) -EE ^~"(g"'^iAi)e-"^'^^"V,07r,(e«®/«).F. 



As for the third summand, one must be more careful, since J^{z) contains zero 
modes, and therefore we have: lim q^'^zJ^(q^"^z) = ±hh(S>t^- This gives us 

m— >+oo 

the following expression for the third term in (3.1): 



-^V.P. X^{q'^Zi/Zj)TTi®TTj{xi(^Xl)J^ 
1=1 m£Z,j 

r 

(3.5) +^Tr($^..$"g-Va;0, 

1=1 

M 

where we define the principal value of a series as V.P. ^ am = lim ^ a^. 
Finally, we obtain the following theorem. 

Theorem 3.1. The correlation function (2.6) satisfies the following system of dif- 
ferential equations: 



(3.6) {k + l)zi—T = - ^fa,h{zi/zj)'ni®'nj{fa.®ea)^ 



aeA+ Jt^i 

a6A+ i/i 
r 

-y^^y^jQ{Zi/Zj)TTi®TTj{xi®Xl)^ 
r Q 

1=1 

where the functions foi,h are defined as follows: 

fa,hit) = Yl ^''{rt)e^''^''\ aeA+ 



f_^,hit) = J2 ^-''iQ'^t)e-^^^-\ aeA+ 



mez 



(3.7) /o(t) = V.P. ^°(0), 

mez 

and the operator Qh equals 

(„mg 
aeA+ ^ 

aeA+ ^ 



+ 2(1 -9-9)'" 



-X 



The functions /a, ft, can be easily written in terms of the theta-functions with 
characteristics. We do not give these expressions here since we are not going to use 
it in this paper. 

Let us look at equations (3.6) more closely. This system looks very much like 
some (twisted) elliptic version of Knizhnik-Zamolodchikov (KZ) system. The only 
term which breaks this analogy is ^ Tii{xi)-^T . Due to this term, this system of 
equations is not closed, since we have no equations which would allow us to find 
the derivatives 

Let us move further. Note that our correlation function T depends not only on 
Zi but also on the modular parameter q. Let us derive the equations for dqT. 
Obviously, 

(3.9) q^T = Ti:i^'- ...^'^q-^e^i-d)) 

oq 
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Now let us substitute into this equation the expression for d given by the Sug- 
awara construction (1.3). This yields: 



d ( 1 \ 

(3.10) {k + l)q-T = E ( E + 2 E ) ' 

where 

am,cc = T>($^ • • • $''g"^e^ : e«[|a| + m£?]/a[-|a| - mg\ :) 

6^,/ = 'n-($i . . . ^""q-^e^ ■ xi\mg\xi\-mg\ :) 

Let us calculate a^.a and h.m,i- Consider first case m > 0. Then, using the 
defining commutation relation for intertwiners (2.1), we get 



= ,l.«l+-^IV($i . . . (.,) . . . $>-^eVa[-|«| - m^]) 

3 

+ TV($i . . . $"e«[|«| + m^]?-^eVa[-|«| - m^]). 
Since ea[/]g-^e'^ = g'e-'^(")g-^e'^ea[/], and [e^f/], = /i^ + //c, we get 



+ . . .$.^^^^.(^,) . . .$-,-^eVa[-|«| - m^]) 

Therefore, denoting (^'^^■^9 q-K<^) = c, we see that for m > 



+ c { k{\a\ + mg) + ) J^. ) . 



"a 



Similar arguments show that 

i 

This gives us the final answer: for m > 0, 
(3.11) 



( - (13^ E(^i/^^)'"'^"^'^^ ® ® e,) + ^ (k{\a\ + mg) + ^) 
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where c = qW\+mg ^-h{ot) ^ 

Similar considerations for m < show that 
(3.12) 



«m,a = i^jY^ |^(%- AO'^'+'^TTi ® ® e«) + ^ [k{\a\ + mg) + 

where c is given by the expression above. 

Now we change the notation: c = . Then the expressions for hm,i are obtained 
in a similar way: 

for m > 0: 

(3.13) hm,l = ^JY^^ Y.^Zj/Zi)'^^7Ti ® 7Tj{xi ® X/) + J^/cW^j ^ 

for m < 0: 

(3.14) = ^ ^^^^^2 ^'^^ihiT^T^i ® ^ji^i ® ^i) + Yzr^^^9^ 
for m = 0: 

(3.15) 6o,/ = Tr{^' . . . $>-^e'^a;^) = ^J^. 

This completes the calculations. We can rewrite the result using the following 
identity: 

(g.)l-l+-^e-M.) _ ..icial/. , (h{^ _ \Ar A 



where z = e^'^^^/^, = e^'^''^/^ and the function is given by 

g27rimrg— 27ria;g27rimC 
]^ g27rimTg— 27ria;\2 ' 

m6Z ^ 

It is easy to check that for < Im ( < Im r we have the following expression for 

mail 



g27rimT g— 2mx 



1 d /^i(7r(a;-C))^'i(0) 



47raa;V ^i(7ra;)^i(7rC) 

This formula defines analytic continuation of ^(x, Q for all C; ^ 7^ ™^ + Recall 
that by definition 

6»i(C|t) =2e^^^/^sin C H ~ e2™^e2'^)(l - e2'^^'^^e-2'^)(l - e^^^^^^) 

n>l 

00 

=2E(-l)'^e'^^^("+^^'sin (2n+l)C. 

n=0 
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Similarly, one can show that 

y ^ ^, = lim ((l>{x, C) + , ■ , )■ 

Expanding the right-hand side of this equation in Laurent series around a; = 0, we 
get that 

^ {qzT' O'l'jO) 1 

Z^^ (1 - g'«9)2 86'i(7rC) 246'; (0) 12' 

Finally, we can formulate the result in the following form: 

Theorem 3.2. The correlation function T satisfies the following differential equa- 
tion: 



(3.16) 

a€A+ ^ i,j \Z7ri g / 

+ 2 E "E^o*^^-?' ~ Ci)7ri <8)7rj(a;i (8)a;/) +vok\ T 
1=1 \ i,j J 

where Zi = e^'^'^'l^, q = e^^^^/s and 

ij27rimT — 27ria; ^ 2 

1 a fe,inix-C))0[{0) 



A7rdx\ ei{7rx)ei{-K() 



^0^^"^"'''^^ 8^1 W 24^i(0) 12 



l^d.l/; Ua,h 2^ I _ (j\a\+mg^-h{a) 1 _ (,|a|+mgg-/i(a 

m>0 ^ m<0 ^ 



^|a|+mgp-/i(a) 



^«''^= E l_q|a|+mgg-fe(a) (|Q'l+^^)+ E 1 _ o|a|+mgg-Ma) ( 1^1 + 
m>0 ^ m<0 ^ 

m>0 ^ m<0 ^ 
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Remark 1. We could as well write Ua^hi ^a,h in terms of theta-functions, but it is 

not necessary: in the next section we show that the terms containing Uot,hi'^a,h can 
be cancelled by a suitable renormalization of J-". 

Remark 2. Note that the function ^ plays a very special role in this 

whole story. The same function also appears later in the integral formulas for JF 
in the case Q — sl2 (see section 5 below); it also appears in the expressions for the 
functions fa,h{z) (cf. Theorem 3.1) for g = 5[2- 

Let us summarize the results obtained so far. We consider ^ as a function of 
zi . . .Zn,q with values in F = Vi ® . . . ® ® ([)«). Then theorems 3.1 and 3.2 
give us differential equations for of the following form: 

d 

^i— = Ai{zi ...Zn, q)T 



(3.18) q—J^ = B{zi...Zn,q)J^, 

where Ai , B are some operators in V which are defined for | | > . . . | z„ | > | qzi \ . 

Note also that we can easily find asymptotics of ^ as g — > 0. Indeed, it is easy 
to show that as g — 0, 



=e '^^^(fc+i) ^ 

where \1/ =< fA,fc, ^^{zi) . . .^^{zn)vx^k > is the "correlation function on the sphere" 
with values in Vi ® . . . ® Ki- For the case when Vi are highest or lowest weight 
representations, which is our main example, this function is very well studied; in 
particular, it is known that it is well defined for l^il > . . . and satisfies trigono- 
metric KZ equation in this region (see [FR] and references therein); there are inte- 
gral formulas for this function (see [SV],[Ch]). So we consider ^' as a well-known 
function. 

Now standard arguments from the theory of ordinary differential equations show 
that the equation (3.18) together with boundary condition (3.19) uniquely deter- 
mine J^; in fact, we only need the equation for q-§^^- 

Remark. Note that we never used the fact that Mx.^k are Verma modules. In fact, 
all the statements of this section are still valid if we consider any modules from the 
category O instead of Verma modules - provided that the intertwiners exist. 

4. Connection with the heat equation 

In this section we will show, following the ideas of Bernard [Ber] that the equation 
(3.16) is closely related to the heat equation on I^r. This is of special interest for 
n = and n — 1. In this cases does not depend on Zi, so (3.16) is the only 
non-trivial equation for J-". 

Let us denote JFq = Tr|Mo ui^~'^^^)- t)e easily calculated explicitly: since 

the character of Mq^^ coincides with the character oiUQ~ , we see that 
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(3.19) 



[] (1 - ^mp^r -Q _ g-/i(a)^|a|+mg)' 

7Ti>0 a,m 

where the second product is taken over aU a e A, m e Z such that \a\ + mg > 0. 
This can be rewritten as foUows: 

(4.2) 

° m>0 a6A+ ym>0 m<0 

Exphcit calculation shows that 
(4.3) 

d 1 ( ^ g|«|+m3g-/^(«) , V- 1 

^ Q;eA+ \m>0 ^ m<0 ^ 

— r > mq—- 

where t'a,/i, 'f^o are defined in Theorem 3.2. 
Similarly, if x e f), 

(4-4) = I xia)uah \ 

where tta^/j, is defined in Theorem 3.2. 

Now let us come back to Theorem 3.2, which we rewrite in the following form 

k r _ 

where we for brevity denoted 
(4.6) 



a€A+ i,j 

r 



2 

where </>(a:^, C); 0o(C) are defined in Theorem 3.2. Comparing (4.5) with (4.4) and 
(4.3), we get the following statement. 
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Proposition 4.1. 

(4.7) 



Proof. The only not immediately obvious step in the proof is dealing with the terms 
containing Ua,h- As for them, 



Comparing it with (4.4), we see that this equals to 

<gradJ^,grad — >=5_; ^ ' 



^ dxi dxi 



□ 

Corollary. 



(4-8) ( - l^,] ^ = 0. 



dq 2 V J^o 



Proof. Rewrite (4.7) in the form 

. . 1 ^ -^f 9 1 ^ \ 1 



dq 2{k + l) V V^o/ fc + 

and substitute there !F = J^q = TtImq ki^~^^'^)- 

So, finally we have the following theorem: 
Theorem 4.1. 



1 



2{k + l] 



^27ri|<a 



A; + 1 



1/90 



i,j \aeA+ 



^0 



/i(a) \a 



27Ti 



9 



T^jifa «) ea) 



(4.9) 



1 A 

+ oy]0o(Ci-Ci)7ri®7rj(a;i®a;i) — ; 

2 «=i / 



where (f){x,(),(f)o{Q are defined in Theorem 3.2. 
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Corollary. Let n = 0, i.e. = 'TT(q "e ). Then / J^q satisfies the heat equation: 

This equation (for untwisted algebra) was first derived by Bernard ([Ber]). 

Let us consider more interesting example n = = TVIm^ ki^^i^)Q~^^'^)- 
this case J-' takes values in some module V over q. From the weight considerations 
it is clear that in fact takes values in the zero- weight subspace V{0). Therefore, 
the equation (4.9) reduces to 

at At 

Note that (f){x, 0) is an elliptic function of a very simple form. Indeed, from the 
expression 



g27rimTg — 27ria; 



X/ Cl _ e27rimTg-277ia;)2 



it is obvious that (t){x, 0) is an even elliptic function with periods 1,t and poles of 
order 2 at the points mr + n, m, n e Z. Finding the leading coefficient at m = 0, 
we see that 



(4.11) 4>{x,0)^-^{p{x) + c), 

for some constant c = c(g) (not to be confused with c of section 3); we won't need 
an explicit expression for it. We will, however, use two properties of c{q): as — > 0, 
c(g) ^ and if g e R then c{q) e R. 
This implies 

Proposition 4.2. Let n = 1. Then the correlation function T satisfies the follow- 
ing equation: 

(4.12) 



f d 1 ^ \ 1 v-^ / /rial h{a)\ \ ^ 



aeA+ 



5. Example: g = sl2,n' = 1- 

Let us consider the simplest possible case g = 5(2, n = 1. Let e, /, /?. be the stan- 
dard basis of sl2- We can identify R ~ by a: t-^ —nixh. Under this identification, 
the Killing form becomes — 27r^a;^, and the Laplace operator: A(, = ~2^'^- 

Let Vfj,, e C be an irreducible highest-weight module over 5(2 with the highest 
weight ^a, where a is the positive root of sl2. In other words, the action of f) on 
the highest weight vector of is given by hv^ = ijlv^. 

Note that the zero- weight subspace Vp(0) is not empty iff /i = 2m, m e For 
this reason, from now on we assume that 1/^ = V2m is a finite-dimensional module 
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of dimension 2m + 1. Then V2m(0) is one-dimensional and fa^a\v2m(o) — ^{m+l). 
Also, it is well-known that in this case the intertwining operator $ : M\^k 
Mx^k ® Vfj. exists and is unique up to a constant factor for generic values of A, k. 
So we can rewrite (4.12) as follows: 



(5.1) f -47ri(/c + 1)^ + ^ J F = m{m+ 1) (p + ^) + c) F, 

where F = J^/J^o is a C- valued function of a; e R, r e 7^, 7i being the upper 
half-plane. 

Note that this equation docs not depend on the highest weight A of the Verma 
module we consider: the only dependence is the boundary conditions (3.19). Note, 
however, that F has the following periodicity property: F{x -|- 1) = e~'^^^F{x); in 
particular, if A e 2Z then F{x) is periodic: F{x + 1) = F{x). 

Note that in the limit A; — > — 1 the equation (5.1) becomes the Lame equation 
(see [WW]): 

(5.2) —F = m{m + 1) [p + ^) + c) F. 

We discuss this relation in the next section. 

To find an explicit expression for F, we use formula (2.5), which allows us to 
rewrite F in terms of untwisted intertwiners $(C) as follows: 

q) = ^> 

where K = 2k,A = X + k. 

This enables us to use well-known integral formulas for ^(C)- For this purpose 
we recall the Wakimoto realization of SI2 ([Wa]), following [BF]. Let us introduce 
the algebra A, generated by 7„, n e Z with the relations: 

[ttn, ttm] = 2nSn+m,0 

(5.3) 

[Pm 7m] — ^n+m,0 

and all the other commutators vanish. 

Next, we can define the module H\, A e C over A as a, module generated by a 
vacuum vector va with the properties: 

a^VA = 0, n > 

.X PnVA = In+lVA = 0, 77, > 

(8.4) ^ 



where k ^ K + 2 ^ 0. 

Define the normal ordering in A by: 
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(5.5) : /?n7m := 

Finally, introduce the free bosonic fields as follows: 





m 


> 




m 


< 




m 


> 




m 


< 



n£2 



(5.6) 

nez 



nez 



Then we have the following proposition, due to Wakimoto [Wa] . 

Proposition, cf. [BF] The following formulas give an action of sh in H^- More- 
over, H\ endowed with this action is isomorphic to Verma module M.k,k for generic 
K.A. 



Jeiz) = Piz) 

(5.7) Jh{z) = -2 : P{z)^iz) : +^/^aiz) 

Jf(z) = - : j^{z)P{z) : +yG^a{z)-f{z)+K-f'{z) 

We also need explicit formulas for the intertwiners in Wakimoto realization. Let 
us introduce the operator e^" : H\ —>■ H^j^2c^ which maps f a to f a+2c^/k ^^'^ 
commutes with all the generators of A except ctQ. Then we can introduce the 
vertex operators X{ji)^ji e f)* by: 



X{ca, z) =: exp [ c / a{z) dz | := exp | c —^z " j exp | c — — 

^ ^ V n<0 / V n>0 ^ 



and the screening operators U (t) : 

U{t) = P{t)X{-^,t) 

Then the intertwining operator ^{z) : Aix^K M.v,K ® can be written as 
follows (we assume that // e Z, so is finite-dimensional, and \ -\- jj, — v — 2m) 



(5.8) ^{z)v = E ( / ^(^«)(-7(^))"f/(^i) • • • U{tm)dt^ ...dt, 

n>0 \ A ^ 



w(8) — 
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where is the lowest-weight vector in and the cycle of integration A C 
C™" \ {ti = 0, z,ti = tj} is chosen so that one can define single- valued branches of 
log ti, log (ti — tj), log {ti — z) on A. General results of the homology theory of local 
systems (see [SV]) show that such cycle is essentially unique. 

In particular, this shows that the zero-mode component ^[0] : M.a,k — M.a,k® 
V^{^) is given by: 

(5.9) $[0] = J X{-^a, z){-^{z)rU{ti) . . . U{t^)dti ...dt^ 

A 

To find F, we must calculate Tt\ma,k ^^[0]q'~^^~'^e~'^'^'*^ . Traces of this type 

can be calculated in the same way as we did in section 3; we borrow the answer 
from [BF, formula 3.14], which can be rewritten as follows: 



\ i=l ^=0 V« 

(5.10) 

ChHj,x[[tl [[ ( 1 xm\[[G{to,ti\^,x) 

i=0 0<i<j<m ^ 1\ \ ) / 

where 

q = e"^" 

^ 2 t^~^ > g27rinx 

G{tQ,ti\T,x) = ^ 2 - o2re-27rix = XI TZ^ 



Applying this formula in our case, we get 



TrI 



//mm \ 
TVIhJ (-7W)"^n/^(^^)n^(-;^'*^)^"''^"'''^"^*M dt^---dtm 

1=1 ^ IV I / / \<i<i<m ^ IV I / / 

Since 
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we finally get 



F{x,q) = z'^i^—^)/'^ 

1 m 
J^,(-I-A)/-.-Q 



- ,_l-A)AfV/»l(T(Co-C,:)|r)^-^""" 



A 



(5-11) n (^^^Wk^) ™' n G{z. U\r,x+\)\dH... dtm 

Since we know that this function in fact does not depend on z (though it is not 
immediately obvious from the formula above), we can put z = 1. In this case we 
can use the following formula: 

^ ' ^ ' ' ' ^ 2 ^i(7ra;)^i(7rC) 

So we have proved the following theorem: 
Theorem 5.1. For any X, k & C, k ^ 0,171 & Z+, the function Fx^i^{x, q) given by 



Fx,.{x,q)^e-^'^^+^'^^ 
J ffje^'^'^^(-i-^)Z(a|T)-2-/« 



A 



^1=1 



(5.13) n '(c. - oir) n g,(.(.+i)(>.(.c.) I ■ ■ 

where 

e^{nC\T) 



q = e""'^ 



satisfies the partial differential equation 



(5.14) (-^^^'^a:^ + -Q^) = + 1) (i + ^) + ^) ^^'^ 
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and the following boundary conditions when q ^ 0, 

(5.15) FxA^, q) ^ C(A, «)e--^+--^V2.^ 
where 

C{X,k) = j |fje-2'^^^*^/'^(sm ttC)"^™/'' 

m \ 

(5.16) n (^^« ^(c^-o))'/'^n 713^2^ 



Remark. The function -F(a:, q) in the theorem above differs from the function F{x, q) 
given by (5.11) by a constant factor. 

Note that in the Umit q ^ the right-hand side of equation (5.14) tends to zero 
and the equation itself tends to the heat equation: 



(5.17) (-27ri«- + ^)/(x,T) = 0, 

and the asymptotics (5.15) of the functions Fa,^ given by Theorem 5.1 form a basis 
in the space of solutions of the heat equation. This allows us to find all the solutions 
of equation (5.14) having good asymptotics as r — > ioo. Let us consider the case 
then T = it, t e M+ and n e iR. In this case equation (5.14) can be considered as 
non-stationary Schrodinger equation: 



(5-18) ^ = J-Hf, 

^ ^ dt 27TK 

with the Hamiltonian 

(5.19) H= — -m{m + l) (p(x + ^) + c(t)) . 

Note that this Hamiltonian is self-adjoint with respect to the norm. 
In the limit q ^ 0, (5.18) tends to 



2 
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Theorem 5.2. Let us fix k & zTR and asssume that t = it,t & W. Let f{x,T) be a 
solution of (5.18) such that 

f{x,r) = fo{x,r){l+g{x,T)) 

where fo{x,T) is a solution of (5.20) and sup |5f(a;, t)| Q as q ^ 0(t +oo). 

Also, let us assume that for some tq fo{x,To) is a rapidly decreasing function ofx. 
Then f{x,T) can be presented in the following form: 

(5.21) I ^^,(,),, 

K 

with the function p(A) given by 

R 

Proof. The fact that Hamiltonian (5.19) is self-adjoint with respect to the L^ norm 
imphes that if f{x, r) is a solution of Schrodinger equation (5.18) then \f{x, t)\l2(_^^ 
does not depend on r. The same applies to /q. Next, /o - as any solution of (5.20) 
- can be written in the form 

fo{x,T) = y e-'^^^^+'^^^^'/2«^(A) 

for some density function p. Letting r ^ and applying inverse Fourier transform, 
we see that 

pW^I J e-'^^fo{x,Q)dx 

R 

Now let us define the function / by 

Then Theorem (5.1) implies that f{x,T) satisfies Schrodinger equation (5.18). 
Moreover, it is easy to check that convergence of F\^k to asymptotics (5.15) is in 
fact uniform in x. Therefore, as q ^ f{x, r) can be written in the form 

f{x,T) = fo{x,T){l+g{x,T)) 

with the same function /o as above and g satisfying sup \g{x,T)\ as q ^ 0. 

xeR 

Comparing this with the assumption of the theorem and taking into account that 
the norms \f{x, r) — f{x, r)|, \fo{x, r)| do not depend on r, we see that /(a;, r) = 
f{x,T), which completes the proof. □ 
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6. Lame functions. 



In this section we use the theorem proved in the previous section to find exphcit 
formulas for Lame functions. The basic idea is to let A; —> —1; obviously, in this limit 

the equation (5.1) becomes Lame equation. Unfortunately, we can not just repeat 
all the arguments of the previous sections for k = —1, since in this case (critical 
level) the whole theory fails; we can not even claim the existence of the intertwiners 
between Verma modules and evaluation representations. Instead, we can use the 
integral formula (5.13) for F and then find the asymptotics when k — > —1. 

We illustrate this general idea on the simplest example m = 1, when all the 
calculations can be done absolutely explicitly. In this case the expression for F 
given at the end of previous section reduces to: 

(6.1) 

^ J \0[m) ^i(7r(x+|))^i(7rC) 

A 

In this case one can take A to be the Pochhammer loop in the t-plane, i.e. the 
element ai o ao o a^^ o ^ of 7ri(C \ {0, 1}, where : [0, 1] — C is the loop going 
around 0, 1 in the anticlockwise/clockwise direction. This cycle is the same as the 
cycle aoo in the ^-plane. 

Now, let us assume that k approaches —1 along the real line from below: k e 
M, /c < —1. Then k < and therefore —2/k>0. In this case we can use the 
saddle-point method for finding asymptotics of the integral (6.1) (see, for example, 
E.T. Copson, "Asymptotic expansions", Cambridge Univ. Press, 1965). 

Proposition. Let f{iy) be defined by the integral 
(6.2) fiiy)^ [ e^^^'U{z)dz, 



A 

where A is a contour in C, w, cj) are analytic functions on A. Let us assume that 
there exists zq & /S. such that: 

1. w'{zq) = ^,w"{zo)^Q 

2. (l>{zo) ^ 

3. For any z E A, z ^ zq we have |e™*^^)| < |e™*^^")|. 
Then f{y) has the following asymptotics as v ^ +00; 



(6.3) f{v) ~ (/)(^o)e^™(^") 



-27r ^ 2 



^^"{zq) 



Moreover, if4>, w also depend analytically on some parameters xi then the asymp- 
totics (6.3) can be differentiated to get asymptotics for -^f[v). 

Therefore, to find the asymptotics of the integral (6.1) we must find critical 
points of e^^^^d\{'KC)., i.e. points Co satisfying 



(6.4) X=-^^<^ 



^i(7rCo) 
24 



Of course, it is hopeless to find explicit expression for (^q as a function of A; 
however, for some special values of A - for example, A G Z - it can be done; we 
return to this later. Therefore, it is better to consider as given and define A by 

(6.4) . 

Let us assume the following: 

1. ^'/(ttCo) 

2. The contour A, described above, can be deformed so that it passes through 
Co precisely twice, and |e^'^''0i(7rC)| attains its absolute maximum on A at Co- 
in this situation we can apply the saddle-point method, which gives us the fol- 
lowing asymptotics for F as ^ — 1: 

(6.5) F{x, g) ~ (1 - e'-'/^){A{q))-^/^y/^ x B{x, q), 



where 



^ ''^ ^i(7r(x+i))^x(7rCo) y 



-7r6';(0| 


r) 


S^log 0i(7rC) 


k=Co 





(6.6) 



The factor 1 — e^^'^/'* in front of the integral appears due to the fact the cycle 
of integration passes through Qq twice. 

Substituting this asymptotics in the equation (5.1), we get: 



^7ri^^^B{x, q) - 27riK—B{x, q) + q^B{x, q) = 2{p{x + -) + c)B{x, q), 
or - letting k — > - 



(6.7) = 2 ( p(a; + -) + c - 27ri^^^ j B{x, q). 

This shows that the function B{x, q), given by (6.6) is a solution of Lame equation 
with m = 1. Therefore, we have proved the following 

Proposition 6.1. Under the assumptions 1 and 2 above, the function 

^i(rr(x+f)) ' 

where A and Co are related by (6.4), is a solution of Lame equation of degree m = 1. 

Let us return to these assumptions. Consider first the case when g e M, Co = 1/2- 
One can easily check that in this case A = 0, and that the function |^i(7rC|T)| 
reaches at C = 1/2 its absolute maximum on the real line. Therefore, in this case 
the assumptions 1 and 2 are satisfied. Therefore, in this case the function f{x) is 
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indeed an asymptotics for the function F{x,q) and therefore is a solution of Lame 
equation. 

In general these assumptions are not true. However, the result is still valid: the 
above defined function f{x) is a solution of Lame equation. To see it, note that the 
assumprions 1,2 are satisfied for q,(o in some neighborhood of qo e R, Co = 1/2. 
Therefore, for all these g, Co the function f{x) defined above is a solution of Lame 
equation. Since / analytically depends on Co it implies that in fact f{x) is a 
solution of Lame equation for all q, Co for which it is well-defined (though in general 
we cannot claim that it is obtained as an asymptotics of the correlation function 
F{x,q)). So we have obtained the following result (which is of course well-know, 
see [WW]): 

Proposition 6.2. For any A, Co satisfying (6.4), the function 



(6.8) fix) 



-nixJli<^-C0 + ^)) 



ei{n{x+^)) ' 
is a solution of Lame equation of degree m — 1. 

It is an instructive exercise to check this by direct computation. 

Let us consider some examples. If A = 0, it is easy to see that the critical point 

is Co = 1/2, and f{x) = G(l, -1\t,x+ f ) = E '^l+TJ = f dn {2Kx), where the 

k' 

periods K,K' of Jacobi elliptic functions are related with q hj q = e~'^^ . More 
generally, if A G 2Z then the critical point is Co = | + ^"T-i and the corresponding 
function f{x) is again dn (2Ka;). 

Similarly, for A = 1 the critical points are Co = ^ , Co = ^ + f , and the function 
f{x) = sn (2Ka;), f{x) = cn {2'Kx); the same functions are obtained for every odd 
integer A. Note that we have two difi'erent solutions corresponding to the same value 
of A; in the best case only one of them is an asymptotics of the correlation function 
on the torus; so the interpretation of these solutions in terms of representation 
theory is still unclear. 

We have shown that A e Z gives us the classical periodic solutions of Lame 
equation of level m = 1: dn (2Kx),sn (2Kx),cn (2Kx) (which are of course well 
known). The general theory (see [WW]) says that there are no other periodic 
solutions. Therefore, the asymptotics of the integral representation of the function 
F for m = 1 gave us all the periodic solutions of Lame equation of degree 1. 

The same technique can be applied for arbitrary m. Though it is not easy to 
see what kind of contour we must take and check the necessary conditions for the 
apllicability of the multi-dimensional saddle point method. Proposition 6.2 suggests 
the following formula: 

Proposition 6.3. For any m G C = Ci • • • Cm £ C'", A G C satisfying 
(6.10) iX + m————-} ^ , ,^ :rTz = for alu = 1 . . .m, 

^ ' ^i(K») ^^1(^(0-0)) 

the function f{x), defined by 



(6.11) fix) = e--'^- n 



^i(7r(x-Ci + i)) 



.^^ ^i(7r(x+f)) 
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is a solution of the Lame equation of degree m. 

This proposition can be proved by a direct calculation, which shows that a 
function of the form (6.11) is a solution of Lame equation if and only if A and Q 
are related by a system of equations (6.10) (see [WW, §23.71]). It is also known 
that almost any solution - not necessarily periodic - can be written in this form. 
The only exception is the case when the constant c in the Lame equation (5.2) is an 
eigenvalue, i.e. when this equation has a doubly-periodic solution. In this case the 
other solution is not periodic and cannot be expressed in terms of theta functions. 
Therefore we see that our approach gave us almost all solutions of Lame equation. 

This general formula has been known for a long time (in [WW] it is given with 
a reference to Hermite's note of 1872), though, of course, it was never obtained as 
an asymptotics of an integral in the way we do. 

Conclusion 

We've considered in detail the case Q = sl2- As for the case of arbitrary g, it 
is known there always exists an analogue of Wakimoto realization (see [FF]) and 
therefore, integral formulas for intertwiners similar to (5.8) (see [ATY]). This in 
turn implies the existence of integral formulas for the traces of the form (2.6). 
Therefore, it is natural to suggest the following conjecture (for the case n = 1, i.e. 
trace with insertion of one intertwiner) : 

Conjecture. If V is a finite- dimensional representation then the correlation func- 
tion T defined by (2.6) has the following asymptotics as k ^ (with all the deriva- 
tives): 

where Ci, C2 do not depend on x, k (hut may depend on g, \). 

This would immediately imply that the function (p/ satisfies the equation 

which is obtained from (6.12) by letting k — > 0. This equation is a natural multidi- 
mensional analogue of the Lame equation. Moreover, if the highest weight A of the 
Verma module is integral then the solutions we get must possess some remarkable 
periodicity properties and should be considered as natural analogues of Lame func- 
tions. Note that in general V^(0) is not one-dimensional, so ^ is a vector function; 
however, for g = s[„ there exist finite-dimensional modules with one-dimensional 
zero-weight component. In this case we obtain the elliptic analogue of Calogero- 
Sutherland operator for the root system A„ (cf. [OP]). This will be discussed in 
detail in our future paper. 

It is an interesting problem to consider the correlation functions in the case where 
V is not a finite-dimensional representation. For example, for = SI2 one can define 
the module = {{xiX2)^ f{xi., X2)\f is a Laurent polynomial in xi, X2 of total degree zero} 
with the action of £[2 given by 

a , d d , d 

e = xi- — , h = xi- X2^ — f = X2^ — . 

(JX2 OX\ 0X2 OXi 
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In this form this module has an obvious generahzation to the case q = sin- 

It is easy to see that W^{0) is one-dimensional and e/|vi/^(o) = Also note 

that if /i e then this module has a finite dimensional submodule of dimension 
2^+1. 

Therefore, the traces of the form (2.6) for the intertwiners ^{z): Mx^k — ^ -^A,fc ® 
Wfj. satisfy the parabolic PDE (5.14) with the constant m = /j, - an arbitrary 
complex number. It is an interesting open question if these correlation functions 
have asymptotics of the form 

K)(/)(g,A,a;). 

Then the function 0(g, A, x) would be a solution of the Lame equation. This would 
provide a method of constructing Lame functions with arbitrary complex value of 
the constant m from representation theory of SI2 
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